Mapping the B,T phase diagram of frustrated metamagnet CuFe02. 
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The magnetic phase diagram of CuFe02 as a function of apphed magnetic field and temperature 
is thoroughly explored and expanded, both for magnetic fields applied parallel and perpendicular to 
the material's c-axis. Pulsed field magnetization measurements extend the typical magnetic stair- 
case of CuFe02 at various temperatures, demonstrating the persistence of the recently discovered 
high field metamagnetic transition up to Tn2 ~ 11 K in both field configurations. An extension 
of the previously introduced phenomenological spin model used to describe the high field magne- 
tization process {Phys. Rev. B, 80, 012406 (2009)) is applied to each of the consecutive low-field 
commensurate spin structures, yielding a semi-quantitative simulation and intuitive description of 
the entire experimental magnetization process in both relevant field directions with a single set of 
parameters. 

PACS numbers: 75.30.Kz,75.10.Hk,75.25.-j 



I. INTRODUCTION 

One of the richest and most fascinating phenomena in 
magnetic systems, geometrical frustration, occurs when 
the specific geometry of an atomic lattice prevents, or 
frustrates, the simultaneous minimization of all mag- 
netic exchange interactions within the system, thereby 
inducing a large magnetic degeneracy. With the primary 
interactions of the magnetic system unable to select a 
unique magnetic ground state, the magnetic behavior 
of frustrated systems is dominated by secondary, often 
weaker interactions, which can vary strongly even across 
closely related materials. Consequently, the field of frus- 
trated magnetism is characterized by its vast richness 
and diversity, exotic magnetic states and low tempera- 
ture physics-l^ilS One of the classic geometries in which 
this phenomenon readily manifests itself is the triangu- 
lar lattice with antiferromagnetic interactions. In ab- 
sence of significant secondary interactions, classical spins 
on a triangular lattice antiferromagnet (TLA) compro- 
mise in their 'desire' to align antiparallel and adopt a 
noncoUinear 120° spin configuration at low temperatu res, 
an underconstrained, highly degenerate ground state.!^ 
The situation can be quite different, however, in sys- 
tems where secondary interactions are significant, such 
as in the stacked delafossite material CuFe02, which 
consists of hexagonal Fe'^"'", 0^~ and Cu"*" layers (space 
group R3m, a = b = 3.03A, c = 17.17A). As the 
Fe'^"'" (3(i^, S — 5/2) ions are the system's only mag- 
netic constituents (Cu+ and 0^~ have filled electronic 
shells) , and their spins interact antiferromagnetically, the 
magnetic system corresponds to an archetypical TLA at 
room temperature (Figurejl^). Strikingly though, in con- 
trast to other delafossite TLAs like LiCr02, AgCr02 and 
CuCr02,'^ CuFe02 adopts a coUinear ground state at 
low temperatures. Based on the electronic configura- 



tion of the Fe^"'" ion (^55/2, L = 0), the antiferromag- 
netic exchange interactions within the system are ex- 
pected to be isotropic, thus yielding a pure Heisenberg 
TLA. The presence of a substantial spin-lattice coupling 
in CuFe02 (the secondary interaction) however, induces 
a low-temperature structural distortion through the 'spin 
Jahn- Teller' effeclP^^, hereby reducing the spin state de- 
generacy in the system. The structural symmetry of the 
CuFe02 lattice is first lowered from the hexagonal RSm 
space group to the monoclinic C2/m space group at Tmi 
w 14 K, to be further reduced to a lower monoclinic sym- 
metry at Tn2 « 11 K.GSHi^l Magnetically, CuFe02 under- 
goes a transition from its paramagnetic (PM) phase to a 
partially disordered, incommensurate (PDIC) magnetic 
phase at Tjvi where a sinusoidally amplitude-modulated 
magnetic structure with a temperatu re dep endent propa- 
gation wave vector (q q 0) is adopted.'SESl Another mag- 
netic phase transition at brings the system into its 
coUinear four-sublattice (4SL) ground state, in which the 
spins align (anti-)parallel to the c-axis, adopting an in- 
pl ane tw o-up two-down order, as illustrated in Figure 
j^j-jjUmi avoid confusion, we will refer to crystallo- 
graphic directions using the hexagonal description de- 
picted in Figure [T] throughout the paper. 

The stabilization of the coUinear 4SL state in CuFe02 
proved to be one of its most puzzling issues. Initially, 
the system was described as a two-dimensional (2D) Ising 
TLA with exchange interactions up to the third nearest- 
neighbors. The first (Ji), second (J2) and third (J3) 
in-plane nearest-neighbor interactions were estimated to 
comp are a s J2/J1 ~ 0.5 and J3/J1 ~ 0.75 in this 
]-^Q jg^isi iTHlo]^ w ith Ji corresponding to approximately 
1.2 meV.'i^'2iJ22| xhere is, however, a priori no physical 
justification for the assumed Ising nature of the mag- 
netic moments. Such an assumption is also inconsistent 
with magnetic susceptibility measurements, which show 
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FIG. 1: (Color online) a) Magnetic structure of CuFe02, space group R3m, a — b = 3.03A, c — 17.17A. Only the magnetic 
Fe"^"^ ions (3ci^, S — 5/2) are depicted, illustrating the magnetic structure of quasi-separate triangular layers. Different 
triangular symmetries as consecutively occurring in CuFe02 are depicted on the right, b) Successively adopted spin 
structures in the various phases of CuFe02, when subjected to an increasing applied magnetic field B \\ c. The higher field 
magnetic phases are proposed on the basis of a recently reported classical spin model (PCS model, see text), c) Analogous 
sequence of consecutively adopted spin arrangements in CuFe02 for the i? _L c configuration. 



highly isotropic behavior above Tjvi in CuFe02 .HSlSSlUlj 
Nonetheless, the magnetic properties below Tjvi are un- 
mistakably strongly anisotropic. The recent discovery of 
the low temperature structural distortion offers an al- 
ternate picture, as it results in a lattice of scalene tri- 
angles in the basal plane (see Fig. [T^), which splits 
the first nearest-neighbor interaction within every tri- 
angle into three unequal exchange interactions, lower- 
ing the energy of the 4SL statcliSHilj Perhaps more im- 
portantly, the distortion has been argued induce a small 
easy axis anisotropy along the c-axis as well, further sta- 
bilizing the collinear ground state.l^ Experimentally, a 
small single-ion anisotropy interaction was estimated by 
fitting a 3D Heisenbcrg Hamiltonian with a single-ion 
anisotropy term to the spin-wave dispersion along the c 



axis below T/vi , wh ich su pports the picture of distortion- 
induced anisotropy.l^i^ As will be confirmed below, the 
combination of this weak magnetic anisotropy and the 
relatively strong spin-phonon couplin g iii CuFe02 can 
explain its observed Ising-like behavior.EEH 

Arguably the most fascinating physical properties arise 
when CuFe02 is subjected to an external magnetic field 
below T]^2 ■ Upon increasing applied magnetic field along 
the c axis (B \\ c), the material has been shown to un- 
dergo a series of magnetic transitions at B^-^ ~ 7 T, 
5^2 ^ 13 T, B^^ ~ 20 T, bI^ ~ 34 T and bI^ ~ 53 
T, b efore ultimately reaching saturation around B^^^^ ~ 
7Q rp| i8 | 22 | 24 | 25 l2g{3j Corresponding magnetic structures 
between the successive transitions (see Fig. [T]d) were 
shown to be a proper helical magnetic order with an 
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incommensurate in-plane wave vector, which also car- 
ries a ferroelectric momentP^IMSl (^ll^ < 5II < j^l^, 
FEIC), and a coUinear five-sublattice (SSL) phase where 
the spins again align (anti-)parallel to the c axis, adopt- 
ing a three-up two-down order {B^^ < i?" < bI^)W^ 
Spin structures at higher fields have not yet been exper- 
imentally determined due to the demanding experimen- 
tal requirements. In a recent work, we have reported 
pulsed field magnetization measurements, revealing the 
retrieval of virtually isotropic magnetic behavior above 
an additional phase transition at JSjlg A corresponding 
anomaly was subsequently observed at somewhat lower 
fields in pulsed-field ultrasonic velocity measurements by 
Quirion et ali^, confirming its proposed magneto-elastic 
nature. On the basis of a phenomenological classical spin 
model (PCS), the spin structures in the high field mag- 
netic phases were suggested to correspond to a coUinear 
three-sublattice (BSL", sjlg < i?" < B^^), an anisotropic 
canted three-sublattice (cBSLH, B^^ < < bI^), and 
an isotropic canted high-field magnetic order (cHF", i?" 

> Bcs), as depicted in Figure [i]d.^^ 

Illustrating the low temperature anisotropy in the ma- 
terial, the magnetism in CuFe02 evolves quite differently 
when it is subjected to a magnetic field perpendicular to 
the c axis {B _L c) below Tjv2, showing only two transi- 
tions up to 40 T, at B^^ ~ 24 T and B^^ ~ 30 TOMmm 
Our recent results also revealed a high field magnetic 
transition for this field configuration, at ~ 51.6 T.^^ 
Apart from the low field 4SL structure, the correspond- 
ing magnetic structures have not yet been experimentally 
determined. Based on the magnetization measurements 
and the aforementioned PCS model, the magnetic struc- 
ture has been proposed to undergo consecutive spin rear- 
rangements from a canted 4SL order (c4SL-'-, with spins 
tilted away from the c-direction) to a coUinear 3SL phase 
(3SL^, with spins in the basal plane) at B^i, to a canted 
3SL order at B^ (cBSL"*-), and finally to the isotropic 
canted high field configuration (cHF-'-) at B^ (see Fig. 

As is clear from above disquisition, the magnetic be- 
havior of CuFe02 as a function of temperature and ap- 
plied magnetic field has proven very rich and has yielded 
unanticipated, fascinating new insights. Following our 
recent results, this work aims to thoroughly map out and 
extend the intricate B,T phase diagrams of CuFe02 up 
to 58 T and Tn2 ~ H K, for both for the B \\ c and 
the i? _L c configuration. Furthermore, by applying the 
recently introduced PCS model to all commensurate sub- 
lattice phases occurring in CuFe02 , an adequate descrip- 
tion of the entire experimental magnetization process in 
both field configurations and an intuitive understanding 
of the magnetic behavior in CuFe02 is provided. 



II. EXPERIMENTAL 

A. Sample preparation 

A high quality, single crystalline rod of CuFe02 was 
synthesized using the Uoating zone technique, following 
the procedure described by Zhao et al.^^ A ^^Fe-enriched 
starting material (^''Fe203, ^''Fe > 95.5%) was used in 
the synthesis, to facilitate nuclear forward scattering ex- 
periments described elsewhere.!^ X-ray Laue diffraction 
was employed to orient the CuFe02 single crystal. Next, 
small cuboid samples (5x1x1 mni"^), with long sides 
oriented parallel (35.9 mg) and perpendicular (42.1 mg) 
to the crystallographic c-axis, respectively, were prepared 
from the single crystal. Further characterization, includ- 
ing ^^Fe Mossbauer spectroscopy, Raman spectroscopy 
and SQUID magnet ometry, also yielded experimental 
data in excellent agreement with literature on CuFe02, 
confirming the high sample quality. The same samples 
were used in all measurements reported here and in a 
previous work.l^ 

B. Instrumentation 

High (pulsed) magnetic field magnetization measure- 
ments, up to a maximum field of 58.3 T were performed 
at the 'Laboratoire National des Champs Magnetiques 
Pulses' in Toulouse, France. The obtained magnetiza- 
tion data were accurately scaled through a least squares 
fit to low field measurements (up to 10 T), performed on a 
well calibrated static (dc) magnetic field setup (using the 
extraction technique) of the 'Institut Neel' in Grenoble, 
France. The accuracy in the scaling procedure was such 
that it introduces an uncertainty of ± 0.3% in all magne- 
tization values determined from the pulsed field experi- 
ments. The temperature dependence of dc magnetic sus- 
ceptibilities of oriented single crystalline CuFe02 cuboids 
was measured in various constant magnetic fields (up to 
7 T) using a Quantum Design MPMS magnetometer. 

III. RESULTS AND DISCUSSION 
A. Magnetization in pulsed magnetic fields 

1. Parallel field configuration (B \\ c) 

Figure [2] depicts the magnetization curves up to 58.3 T 
for various temperatures below T/vi, where the applied 
magnetic field B is parallel to the c-axis (S"). As i?" 
increases, several successive metamagnetic steps are ob- 
served, hiexcellfint agreement with previously reported 
j.ggyl^gpl22l2ll25l27H3il At 1.5 K, the system shows mag- 
netic phase transitions at sj.'^ ~ 7.2 T (4SL to FEIC 
phase transition), bI^ ~ 13.0 T (FEIC ^ 5SL), b]^ ~ 
19.7T (5SL^ 3SLII), bI^ ~ 32.4 T (3SLII cBSL") and 
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B^'g - 53.3 T (cSSLII cHF"). The three transitions at 
lowest critical fields B^^, and B^^ are all accompa- 
nied by large magnetization steps and exhibit significant 
hysteresis (B^^^ = 7.27 T, B^^^ = 7.15 T, bI^^ = 13.44 

T, bI^^ = 12.51, 5^3^ = 20.32 T and B^,^ = 19.08 T 
at 1.5 K), indicating their first order nature. In con- 
trast, at the fourth magnetic transition (B^^ ~ 32.4 T), 
the Af (S'l )-curvc shows only a change in slope, suggest- 
ing this transition is of second order, whi ch is c onsistent 
with synchrotron x-ray diffraction results. '2^1221 xhe high 

field transition at is again of first order nature, as 

illustrated by its hysteresis: B^^^^ = 53.78 T and B^,^,^ 
= 52.88 T at 1.5 K. The existence of this high field 
transition was recently confirmed in ultrasonic velocity 
measurement^^ and can also be seen in previous magne- 
tization data recorded by Ajiro et al, who measured the 
magnetization of a powder sample of CuFe02 at 8 K in 
a single turn coil measurement up to 100 T.l^ Though it 
is obscured in their A/,_B-curve, presumably due to the 
polycrystalline nature of the sample, a clear feature can 
be seen around ^ 52 T in the corresponding (dM/dB) 
vs. B graph. 

In the 4SL phase, the magnetization is close to zero, 
as expected for the two-up two-down structure (Fig. [T|d, 
4SL). In the FEI C phase, M increases linearly with ijH 
as observed befor^^ZHmiElIST]^ signaling a continuous re- 
orientation of the spin system in the spiral phase (Fig. 
[T|d, FEIC). In the 5SL phase, M is almost constant, at a 
value approximately equal to one- fifth of the 5 fj,B /Fe^'^ 
saturation value, in good agreement with the three-up 

two-down structure (Fig. Ip, 5SL). Between bI-^ and 

II I 
5^4, M is again almost independent of , having a value 

close to l/3rd of the saturation- value, while between B^^ 
and i?),'^ the magnetization again increases linearly with 
-B", indicating another continuous reorientation of the 
spin system. Based on these observations and the PCS 
model, these phases have been proposed to correspond 
to a coUinear three-sublattice (Fig. ^p, 3SLII, two-up 
one down) and a canted three-sublattice phase (Fig. [T|d, 

c3SLII ), respectively.^'^ At sjlg, the system undergoes 
another first order transition, where the magnetization 
exhibits an abrupt jump. Above B^j.^, the magnetization 
shows a steady linear increase up to the highest field mea- 
sured, 58.27 T. At this point M has taken a value of 3.54 
Hs/Fe^^ (at 1.5 K), close to the ~ 3.7 ^B/Fe3+ value 
for the powder sample measured at 8 K by Ajiro et alP^- 
As the system has regained isotropic behavior above this 
transition, the spin structure in this cHF" phase has been 
proposed to be isotropic, such as e.g. the canted 120° 
configuration depicted in Fig. [T]3, where the projection of 
the spins in the basal plane retains the typical 120° con- 
figuration while the out of plane spin-components grow 
with B". In their recent paper, Quirion et al. proposed 
a similar, though slightly incommensurate 120°-like spin 
structure based on Landau free energy considerations.!^ 



As the temperature increases, the general features of 
the Af ,i?-curve survive, though magnetic steps are broad- 
ened over an increasingly wide field range, hysteresis 
widths are reduced and plateau phases acquire increasing 
slopes. As the temperature approaches T/V2j the charac- 
teristic staircase features of the magnetization smooth 
out and M increases (quasi-)linearly with B, deviating 
from this behavior only at high magnetic fields, close to 
saturation. The fact that this appears to occur already 
just below T/V2 is ascribed to a slight offset of the cor- 
responding temperature sensor at these temperatures, as 
transition temperatures measured in susceptibility exper- 
iments on the same sample (see below) are in accordance 
with literature values. A striking feature is the tempera- 
ture dependence of the various magnetic transitions (in- 
dicated by the dashed arrows in Figure[2]). Figure|3]shows 
the relative variation of the corresponding critical mag- 
netic fields with temperature. With the exception of the 
lowest field-induced transition, all (first order) transitions 
exhibiting hysteresis show identical behavior; a continu- 
ous decrease of the corresponding critical field {b\2, 
and Bllg, respectively) with increasing temperature. In 
contrast, the critical field of the second order transition 
(^c4) proves rather temperature independent, once more 
indicating its different nature. 



2. Perpendicular field configuration (B A- c) 

Figure [4] shows the magnetization process up to 58.3 
T for various temperatures below T/vi, for the perpen- 
dicular configuration (i? _L c). As for the parallel config- 
uration, the magnetization cu rves are in ex cellent agree- 
ment with earlier observations! "'^^*^^!^^*^^ ! ^^ * With increas- 
ing B-L, the magnetization shows a steady linear increase 
up to Bjr^ (~ 24.8 T at 1.5 K), suggesting a slight continu- 
ous canting of the 4SL spins from the c direction, toward 
the basal (a, 6) plane (c4SL, Fig. [ij;). Indeed, neutron 
diffraction data have confirmed the stability of this c4SL 
magnetic structure up to at least 14.5 T;-^ At Bjr^, the 
system undergoes a first order magnetic transition to a 
plateau state, which shows significant hysteresis (at 1.5 
K, B^^^ = 25.40 T and B^^^ = 24.27 T). The magnetiza- 
tion in this plateau state is rather independent of B^ at 
an average value of ~ 1.53 /is/Fe^+, close to l/3rd of the 
saturation value, implying a three-sublattice state with 
spins in the basal plane, directed along B^ (Fig. [l];, 
3SL-'-). This spin configuration was recently confirmed 
using numerical minimization of the PCS model. ■^^ Af- 
ter undergoing a second order phase transition at ~ 
30.0 T (at 1.5 K), M once again increases (quasi-)linearly 
with B-^ ^ which in turn implies a continuous reorienta- 
tion of the moments away from collinearity. Due to the 
nonzero easy axis anisotropyat these fields, the slope in 
this canted 3SL phase (Fig. lip, c3SL^) differs from that 
in the same field interval for the parallel configuration. 
At Bjr^ ~ 51.6 T (1.5 K), another magnetic transition is 
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B ' (T) 

FIG. 2: (Color online) Magnetization measurements in pulsed magnetic fields at various temperatures. The magnetic field is 
applied in the direction parallel to the c-axis. The various curves are offset by consecutive multiples of 0.5 fj-s/Fe^^ with 
increasing temperature for clarity. Thick and thin lines represent sample magnetization in increasing and subsequently 
decreasing magnetic field, respectively. Dashed arrows indicate the temperature dependence of the various magnetic 

transitions (See also Fig. |3|. 
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FIG. 3: (Color online) Temperature dependence of the critical fields corresponding to the various magnetic transitions, for 
both configurations {B || c and _B _L c). Critical field values (midpoints of hysteresis loop for first order transitions) are 

normalized by their corresponding 1.5 K values. 



observed, similar to that at B^^ in the parallel configu- 
ration. As in that configuration, the high field transition 
here consists of a first order metamagnetic step, which 
exhibits hysteresis (at 1.5 K, B^^^ = 52.02 T and = 
51.18 T). At 1.5 K, M jumps to ~ 3.1 ^lB|^e^"^ at B^^^, 
after which it resumes a steady increase, in line with a 



noncollinear spin arrangement. In fact, based on the ab- 
sence of anisotropy in this canted HF phase (cHF-'-), the 
PCS model predicts a spin structure analogous to that 
for i? II c, as sketched in Fig. [l|:. The fact that the addi- 
tional transition occurs at slightly lower critical field in 
the perpendicular configuration {Bjr^ ~ 51.6 T vs. b\c^ 



6 



~ 53.3 T) explains the broadness of the feature around 
'--^ 52 T in the aforementioned {dM / dB) ,B-cmve of the 
polycrystalhne sample of Ajiro et al^ 

With increasing temperature, the general features of 
the M,_B-curve remain intact, although the plateau phase 
acquires an increasing slope. Furthermore, as for the par- 
allel case, the transition features are smoothed out upon 
approaching T/v2, and deviation from this behavior only 
occurs upon approaching saturation. Again, the appar- 
ent small temperature mismatch with respect to suscep- 
tibility measurements (below) is attributed to a slight 
offset of the temperature sensor at temperatures close to 
Tn2- Also for B ± c, the temperature dependence of the 
various critical fields correlates to the nature of the corre- 
sponding transitions (See Fig. |3|; first order transitions 
(at B^i and B^) exhibit the same relative decrease with 

temperature as i?|.2, bI^ and B^^, while the second or- 
der transition (at B^) shows a much weaker temperature 
dependence. 



3. Progressive symmetry increase 

The strong coupling between spin and lattice degrees 
of freedom is a key ingredient in the description of the 
m agnet ization process of CuFe02. Recently, Terada et 
oi;j 25 | 30 | gjiowed the strong correlation between the lat- 
tice parameters and the magnetization in applied field 
in both configurations. For B \\ c, coinciding with the 

metamagnetic steps at bIj_, BI2, and sjlg, the lattice un- 
dergoes corresponding discontinuous contractions along 
the [110] direction, while changes in the [110] direction 
are much smaller. In addition, the lattice has been shown 

to increase its symmetry at B^^2^ where the scalene trian- 
gle distortion is partially relieved, resulting in a lattice 
of isosceles triangles (Fig. [l^).'^ The lattice parameter 
along [110] mirrors the behavior of the magnetization in 
applied field; within the collinear phases it remains prac- 
tically constant and in the noncoUinear phases the lattice 
continuously contracts with increasing field (and magne- 
tization) . These striking observations can be rationalized 
as follows: in zero field, the spin-lattice coupling induces 
the scalene triangle distortion and a magnetic easy axis 
along the c direction, thereby reducing the magnetic en- 
ergy at the expense of elastic energy. As i?" increases 
however, a growing tendency for parallel spin alignment 
in the field direction develops, thereby successively re- 
ducing the degree of magnetic frustration (the driving 
force for the distortion). Thus, as the gain in magnetic 
exchange energy is successively reduced with , the sys- 
tem rebalances the magnetic and elastic energies asso- 
ciated with the lattice distortion along with every spin 
rearrangement. As a result, the system exhibits a pro- 
gressive lattice contraction along [110], which mirrors the 
changes in magnetization. 

Since the induced magnetic anisotropy in the mate- 
rial is also directly coupled to the lattice distortion. 



one may expect the strength of the induced single-ion 
anisotropy to diminish accordingly with M, undergoing 
steps across first order transitions and continuously de- 
creasing in (giiasz-) linear phases. Indeed, as shown in 
Figure 2 of our recent papei^, which shows the M, in- 
curves for both the parallel and perpendicular configura- 
tion at 1.5 K, the system's response to an applied field be- 
comes more and more isotropic as B increases. Moreover, 
above both high field transitions, recently confirmed to 
be magneto-elastic in nature, the system was even found 
to behave almost completely isotropic, consistent with 
a vanishing easy axis anisotropy and the retrieval of an 
undistorted equilateral triangular lattice. 



B. Magnetic susceptibility in constant fields 

In order to supplement the magnetic phase diagrams of 
CuFe02 and to further elucidate its magnetic behavior, 
the temperature dependence of dc magnetic susceptibili- 
ties in various constant magnetic fields was measured in 
both field configurations. Panels a)-c) of Figure [5] com- 
pare the low temperature magnetic susceptibility curves 
for the two field orientations in applied fields of 0.01, 4 
and 7 T, respectively. Consistent with previous measure- 
ments, both Xm II ^) ^^'i Xm -L c) show a diffuse 
maximum at Tjvi ~ 13.5 K and a subsequent abrupt drop 
at Tiy2 — 11-2 K upon decreasing temperature.'i^^SSMM] 
Above Tn2 the susceptibility is isotropic, for all applied 
fields measured. As expected for an ordered antiferro- 
magnet, x^m approaches zero with decreasing tempera- 
ture below T/V2j while Xm remains almost constant af- 
ter the initial drop at T/v2. The field dependence of the 
magnetic susceptibility, visualized in panels d) and e) 
for the parallel and perpendicular configuration, respec- 
tively, shows the invariance of T/vi with applied field for 
both configurations. Though relatively field independent 
for the perpendicular configuration, Tn2 shifts to lower 
temperatures as the applied magnetic field approaches 
B^.^ (~ 7.2 T) in the parallel case. This difference can 
be regarded as a consequence of the lower susceptibility 
in the ordered phase for B \\ c, which is unfavorable to- 
ward the Zeeman interaction, which becomes increasingly 
strong with B. Thus, with increasing i?" the magnetic 
ordering transition at tJi^^ is shifted to lower tempera- 
ture. For the perpendicular case, the susceptibility drop 
across is only marginal, ergo the corresponding tem- 
perature down-shift is far less pronounced. 

As is clear from panels d) and e), the transition at rj!|2 
and 7 T acquires a double feature, indicating the process 
becomes two-stepped. This points toward the presence 
of an intermediate phase between the two steps. Based 
on the constructed phase diagram presented below (Fig- 
ure [6| , this intermediate phase is identified as the helical 
FEIC phase, as the phase boundaries of both the 4SL 
and FEIC phases bend toward the Tjl^^i^) line at these 
temperatures. 
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-F/G. J^: (Color online) Magnetization measurements in pulsed magnetic fields at various temperatures. Here, the applied 
magnetic field is perpendicular to the c direction of the crystal. The various curves are offset by consecutive multiples of 0.5 
/iB /Fe^^ with increasing temperature for clarity. Thick and thin lines represent sample magnetization in increasing and 
subsequently decreasing magnetic field, respectively. Dashed arrows indicate the temperature dependence of the various 

metamagnetic transitions. 
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FIG. 5: (Color online) Magnetic susceptibility versus temperature in various constant magnetic fields for both the parallel 
and the perpendicular configuration. Panels a)-c) compare the magnetic susceptibility for -B |! c and B _L c as a function of 
temperature at selected fields up to 7 T. Vertical dotted/dashed lines indicate magnetic transition temperatures, as deducted 
from corresponding dxM /dT data. The zoom-ins on the data in panels i) and j) depict the field dependence of the magnetic 

susceptibility in parallel and perpendicular configurations, respectively. 



C. Phase diagrams 

With the phase transition data obtained above in hand, 
the experimental phase diagram of CuFe02 as a func- 
tion of applied field and temperature can be assembled 
for both field configurations. Magnetic transition fields 



and temperatures are defined through the position (cen- 
ter) of corresponding anomalies in the derivatives of the 
pulsed field magnetization (dM/dB) and susceptibility 
(dxAi/dT) curves, respectively. 

Figure [6] shows the B,T phase diagram for CuFe02 
that can be constructed based on aforementioned 
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FIG. 6: (Color online) Top panel: Differential magnetization curves for B \\ c. Colored solid and dotted lines depict 
(normalized) d A'l /dB -cnrves as measured in increasing and subsequently decreasing magnetic field, respectively. 
dM / dB-curves measured in decreasing magnetic field are inverted for clarity. Additionally, each curve has an offset equal to 
its corresponding temperature . Lower panel: B,T phase diagram of CuFe02 for the case where B is parallel to the c-axis. 

Large, black squares depict magnetic transitions as observed in this work and smaller gray symbols indicate previously 
reported transitions. Open and closed symbols represent transitions observed in increasing and decreasing {B or T) ramps, 
respectively. Solid lines are correspond to proposed phase boundaries. Region / corresponds to an observed intermediate state 

(see text). 



exp eriments and othe r, currently available literature 
^g^^^mmmsmmMSm^ for the parallel configuration 
{B II c). The diagram features all the previously con- 
firmed phases; the zero field PM, PDIC and 4SL phases 
and the consecutive FEIC ^ SSL ^ SSL" ^ cSSL" ^ 
cHF'I phase cascade upon increasing field below Tn2- 
Worth noting is the fact that the transition from the 
5SL to the SSL phase (at B^j.^) appears to split up into a 
tvi^o-step transition with temperature, implying an inter- 
mediate spin state /. At temperatures approaching T7V2, 
the magnetization of the system in the corresponding 
field region deviates continuously from the 3SL plateau 
value (see the 7 K line in Figure [2]), suggesting that here 
(some) spins are canting away from collinearity, before 
the abrupt rearrangement to the SSL spin configuration. 



These double transition features were observed before in 
steady state magnetic field measurements up to 23 li^, 
which indicates that this behavior reflects the inherent 
reduction of the magnetic anisotropy with applied mag- 
netic field in CuFe02. 

Figure [7] shows the analogous B,T phase diagram 
resulting from above experiments and earlier reported 
data^^^^ for the case where i? _L c. The diagram in- 
cludes the zero field PM, PDIC and 4SL phases and the 
field-induced phase cascade for this field configuration: 
c4SL ^ 3SL-L ^ C3SL-L ^ cHF-^. 
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FIG. 7: (Color online) Top panel: Differential magnetization curves for B _L c. Colored solid and dotted lines depict 
(normalized) dM /dB-cvLTves as measured in increasing and subsequently decreasing magnetic field, respectively. 
dM/dB-curves measured in decreasing magnetic field are inverted for clarity. Additionally, each curve has an offset equal to 
its corresponding temperature for clarity. Bottom panel: B,T phase diagram of CuFe02 for the _B _L c case. Large, black 
squares depict magnetic transitions as observed in this work and smaller gray symbols indicate previously reported 
transitions. Open and closed symbols represent transitions observed in increasing and decreasing {B or T) ramps, 
respectively. Solid lines correspond to proposed phase boundaries. 




FIG. 8: (Color online) Effective magnetic unit cells in the consecutive a) four-, b) five- and c) three-sublattice phases of 
CuFe02. In all three structures the dotted black lines outline the magnetic unit cell when a single layer is considered (i.e. 
when Jz = 0), while the grey shaded volume indicates the effective two-layer magnetic unit cell considered when interlayer 
interactions are taken into account. Thus, an energetically optimal ABAB-type stacking is assumed for each sublattice. 



D. Classical spin model 

In order to study the magnetization process in CuFe02 
further, we resort to the PCS model. This phenomeno- 
logical model includes the primary magnetic interactions 
of the system; along with the basic magnetic exchange 



and Zeeman interaction terms, the strong spin-phonon 
coupling and the magnetic isotropy in CuFe02 arc in- 
cluded. The incorporation of the latter two seems key 
to capture the Ising-like behavior of the system, as was 
recently shown.l^I^Zl To determine the effect of interlayer 
exchange interactions on the system, these are included 
in the model separately later. 
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Spin- lattice interactions are typically incorporated into 
the Hamiltonian throug h the d istance dependence of the 
exchange coupling J(r) pifilll] Ergo, in general for a sys- 
tem with isotropic exchange interactions the effective 
Hamiltonian becomes: 

^e//. = 5^ • ^^(l - au^j) + HdefX{u^}), (1) 

where the Ui are the displacement vectors, the Uij 
{=(ui — Uj)-rij /\rij\) are the corresponding relative 
changes in bond length between sites i and j, a is the 
spin-lattice constant (to first approximation equal to 
dJ/dr) and Hdef. corresponds to the deformation energy 
cost associated with the atom displacements Ui, which is 
thus dependent on the phonon model still to be chosen. 
Taking the simple bond-phonon (BP) model here, which 
treats the bond lengths Uij as independent variables, the 
presence of spin-phonon coupling effectively introduces 
an additional biquadratic spin interaction of stre ngth 6 J, 
where b = a^J/k (third term in eq. [2|.EnilIl Further- 
more, since neighboring bond lengths independent 
here, the biquadratic term is restricted to nearest neigh- 
bor couplings only. Due to the quadratic nature of the 
term, either parallel or antiparallel spin configurations 
are favorable, which explains the tendency of spin-lattice 
coupling to stabilize coUinear spin states. 

Thus, the general spin Hamiltonian (containing only 
magnetic contributions) for CuFe02 within the PCS 
modeP^ can now be constructed: 

-Y^bJ.,iS,-S,f-DiB)J2sl, (2) 

where B is the applied magnetic field, Jij is the exchange 
interaction between sites i and j, b is the biquadratic 
coupling constant and D is the magnetic anisotropy con- 
stant, which is field dependent due to its strong coupling 
to the lattice distortion. The Zeeman and anisotropy 
terms sum over all sites i, the biquadratic term couples 
only nearest neighbor spin pairs i and j, and the exchange 
term includes all spin pair interactions in the system. 

In a previous work, we analyzed the behavior of this 
spin Hamiltonian (eq. [2]) when applied to the magnetic 
unit cell of the three sublattice structure, thereby focus- 
ing on the high field magnetic phases of CuFe02. Here, 
we compare the field dependence of all consecutive com- 
mensurate phases, based on the same spin Hamiltonian 
and the previously extracted parameters. Thus, we eval- 
uate the corresponding spin Hamiltonians for the mag- 
netic unit cells of the four-, five- and three-sublattice 
structures on a single triangular sheet; the correspond- 
ing unit cells are sketched in Figure [8j Considering the 
spins as classical, justified by the large S = 5/2 value, 
we write Si = BiS (where e is a unit vector), and Ji, J2 
and J3 for the first, second and third nearest neighbor 
exchange interactions, respectively. The respective spin 
Hamiltonians are then found to be: 



H4SL — 

i i 

+ 2JiS^{pi2 +P13 +Pi4 +P23 +P34) 
+ 2J2S^ (2 + P12 + pii + P2Z + Pm) 

-h4J3S'2(l+pi3+p24) 

- 2AS^{pi2 + p?3 + pii + pI^ + pli + , (3) 
H<:,SL = 

^2^iBSB - Y^e,- D{B)S^Y.^h 

i i 
+2JiS'^{pi2 +P15 +P23 +P34 +P45) 
+ JlS'^ (pi3 + pu + P2i + P25 + P35 ) 
+ 128"^ {5 + 2pi3 + 2pi4 + 2p24 + 2p25 + 2p35) 
-|-2J3S'^(pi3 +_pi4 +P24 +P25 + P35) 
+ J3S'^{pi2 +P15 + P23 + P34 + P45) 

-2AS\pl2+pl5 +P23 +P34 +P45) 

-^^^b?3+P?4+P24+P25+P35), (4) 
HsSL — 

^2tiBSB • ^ e, - D{B)S^ ^ e,?, 

i i 
-F3(Ji + J3)S\pi2 +P13+P23) 

^3AS\pl2+pl3+Pl3): (5) 

where g is taken as 2 and spin-spin couplings are writ- 
ten as pij (= ei ■ Bj). The exchange constants are taken 
as equal along the different in-plane crystallographic di- 
rections, their field- dependence being in the spin-phonon 
term. The spin-phonon parameter is defined as ^ = bJi, 
which corresponds to G/3 in our previous work^^, though 
with a rescaled dimensionless biquadratic coupling b of 
~ 0.0098 (here, third nearest neighbour interactions are 
taken into account in the estimation of b). 

To test the PCS model, we performed numerical mini- 
mization of equations[3][5]as a function of the independent 
spin vectors (ci, 62, 63, 64 and 65) at a given field B, 
using previously extracted parameters. In order to incor- 
porate its field dependence, which is a priori unknown, 
D is approximated to be proportional to {Msat. — M{B)) 
here. In other words, D is assumed to mirror the field- 
dependence of Af , undergoing stepwise reductions at first 
order transitions and vanishing as the system approaches 
saturation; see Figures [lO] and [TT] The previously esti- 
mated D for the collinear 3SL phases (SSL" and 3SL^) 
was ~ 0.021 meV, making it ~ 0.031 meV and ~ 0.025 
meV in the collinear 4SL and 5SL phases, respectively. 
Taking exchange couphngs as Ji ~ 0.259 meV, J2 ~ 0.102 
meV and J3 ~ 0.181 meV, and the spin-phonon param- 
eter A as ~ 0.00247 meV, the parallel-field dependence 
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{B II c) of the resulting energy per spin for each of the 
commensurate sublattice phases is as shown in Figure [9^. 

Upon examination of the different energy curves, one 
finds that the PCS model with these parameters yields a 
cascade of expected magnetic transitions that is consis- 
tent with experiment. The 4SL collinear four-sublattice 
state is stable with respect to the 5SL structure up to 
~ 9.4 T. From there on, the 5SL state is the most en- 
ergetically favorable, up to the critical field bI^. Above 
B^^, the colhnear SSL^ state becomes stable, undergo- 
ing a transition to the cSSL" structure only around sjl^ 
~ 32.4 T. Experimentally, the multifcrroic spiral FEIC 
phase was found as an intermediate phase, between _b||j^ 

~ 7.2 T and — 13-0 T. As this phase is incommen- 
surate, however, it is not feasible to describe it using 
the PCS model applied to a limited-size unit cell here. 
Recently though, such complex incommensurate ground 
state structures were found in zero fiel d for far larger 
unit cells using Monte Carlo simulations.l^^Illl Based on 
the experimental data, the energy per spin of the field- 
induced FEIC phase in CuFe02 is expected to have a 
field-dependence as indicated by the green dotted line in 
Figure [9^, making it the adopted spin structure between 
Bl^ndB^,. 

The corresponding magnetization curve for B || c, de- 
picted in Figure [TO^, shows a good agreement with the 
experimental result (Figure [2) 1.5 K curve). The non- 
directional spin-phonon interaction (A), which favors 
collinear spin states, combines with the directional ap- 
plied field i?" and the easy axis anisotropy D{B) to suc- 
cessively stabilize the consecutive magnetization plateaus 
of the collinear phases. At high fields (above b\^), the in- 
creasingly dominant Zeeman term and the progressively 
reduced anisotropy result in a gradual spin canting in the 
system. 

In analogous fashion, one can calculate the energy per 
spin for the commensurate sublattice phases in case of 
a field applied perpendicular to the c-axis using equa- 
tions [3] and [5j respectively. Using the same parameters 
as used for the B \\ c case, one obtains an energy scheme 
as depicted in Figure [TT^ . The c4SL state is the most 
energetically favored up to -B^ ~ 24.8 T, above which a 
three-sublattice is the most stable, with the spins adopt- 
ing consecutive 3SL-'- and 3SL-'^ structures as the applied 
field increases. The inset of Fig. [TT^ shows the corre- 
sponding simulated magnetization curve for i? _L c, as 
well as the corresponding assumed value of the magnetic 
easy axis anisotropy in the various magnetic phases. The 
obtained magnetization process is again in good agree- 
ment with the experimental curve (Figure |4) 1.5 K line). 
As opposed to the B \\c case, the directional anisotropy is 
orthogonal to the field direction here, resulting in a much 
smaller plateau width. Thus, the PCS spin Hamiltonian 
(eq. [2]) also provides an adequate description of the low 
field part of the magnetization process in CuFe02, for 
both field configurations, using the same parameters that 
were previously used for describing the high field part. 



We emphasize the fact that the spin Hamiltonian pa- 
rameters used were determined through direct compar- 
ison with experimentally observed features. The easy 
axis anisotropy D (only a scaling parameter as D(B) oc 
{Msat. — ^{B))) and spin-phonon coupling A were de- 
termined through the simulation of the high-field mag- 
netization process, which also set the value for the sum 
of Ji and Js.^'' With these preset restrictions, J2 and J3 
were set such that: i. the simulated 5SL to 3SLII transi- 
tion field for i? II c corresponds to the experimental value 
{b\.^), and a. the simulated c4SL to 3SL''" transition field 
corresponds to the experimental i?^ value. The result- 
ing exchange parameters compare as J2/J1 — 0.39 and 

J3/J1 ~ 0.70, r atios which are close to those previously 
estimated.l^ilSei 



E. Interlayer exchange interaction 

Recent efforts suggest a magnetic exchange interaction 
between the Fe-layers to be an additional important as- 
pect of the CuFe02 system. Inelastic neutron scatter- 
ing work shows indicative spin- wave dispersion along the 
hexagonal axis, sign aling the interlayer interaction to be 
significantl ^^ l ^^ l ^^ *^ This is corroborated by the observa- 
tion of finite dispersion of calculated electronic bands.l^ 
Thus, here wc incorporate the interlayer exchange into 
the PCS model to determine its influence on the modeled 
magnetization process. As the interplane interaction is 
estimated to be small compared to the in-plane exchange, 
we take a perturbative approach, taking only nearest 
neighbor interactions (Jz>0). With each spin having 
nearest neighbor interlayer couplings to three consecutive 
sublattice sites in the adjacent layers, all types of stack- 
ing of three-sublattice layers are energetically equivalent, 
while four- and five-sublattice layers have specific opti- 
mum stacking sequences (those depicted in figure [8| .'^SEIl 
Assuming this optimal stacking of consecutive layers, the 
effective magnetic unit cells of the four-, five- and three- 
sublattice structures now contain two triangular sheets 
each, with the additional interlayer interactions amount- 
ing to: 



Hz,4SL 
Hz,5SL 



2JzS'^{pi2 + Pl3 + Pl4 + P23 + P24 + P34I6) 
2JzS^{pi3 +pii +P2i +P25 +P35) 



+ JzS^{pi2 +P15 +P23 +P3i +P45) 



Hz,3SL = JzS^{3 + 2pi2 + 2pi3 + 2p23) 



(7) 
(8) 



per layer. For the collinear sublattice structures these 
terms add up to —J^S^, ~JzS^ and +JzS^/3 per spin 
for the four-, five- and three-sublattice structures, respec- 
tively. With the inclusion of these terms, equations |3][5] 
were once again numerically minimized to determine op- 
timum spin directions in an increasing field; the resulting 
energies of the different sublattices and the correspond- 
ing simulated magnetization curves are depicted in Fig- 
ures |9]d and 10 d for B \\ c and in Figure [Tl| 3 for B _L 
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FIG. 9: (Color online) Calculated minimum energy per spin for each of the commensurate sublattice phases {B \\ c) as given 
by numerical minimization of equations 3|5 without (a)) and with interlayer exchange interactions (b)). Solid black, red and 

blue lines correspond to the four-, five- and three-sublattice structures, respectively. The dotted green line represents the 
expected energy per spin of the incommensurate spiral (FEIC) phase. Dashed vertical lines indicate experimental transition 
fields. Insets show a zoom-in on the 5-15 T region, which features two magnetic phase transitions. 
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FIG. 10: (Color online) Field dependence [B \\ c) of the magnetization and easy axis anisotropy in CuFe02 for a) = 
and b) Jz = 0.25Ji. Solid blue line: simulated magnetization process of CuFe02 at low magnetic fields (see text). The blue 
dotted line in the FEIC phase corresponds to the experimental data at 1.5 K (increasing field). Solid red line: corresponding 
assumed values of the magnetic easy axis anisotropy in the simulation of the various magnetic phases; D is approximated to 

be proportional to (Msat. — M(B)). 



c. As is clear from these graphs and their comparison to 
the case where Jz = 0, the experimental magnetization 
process is equally well simulated upon incorporation of 
interlayer interactions. Keeping A and D{B) at the same 
value, the incorporation of J^, which was fixed at 0.25Ji 
(a represe ntativ e value based on inelastic neutron scat- 
tering datcPillSl)^ results in adapted extracted exchange 
couplings of Ji ~ 0.215 meV (making ~ 0.054 meV), 
J2 ~ 0.068 meV and J3 ~ 0.195 meV. As before, these 
parameters were determined through direct comparison 
with observed experimental features of the magnetization 
process. Though the introduction of an additional anti- 
ferromagnetic interaction in the model generally tends to 
decrease the extracted parameters, J3 is in fact increased 
here to counter the relative destabilization of the three- 
sublattice structure. Summarizing, incorporation of in- 
terlayer interactions into the PCS model yields an equally 
adequate description of the experimental magnetization 
process of CuFe02, with slightly modified exchange pa- 



rameters. 

At this point, it is worth pointing out the limitations 
of the PCS model presented here. As our calculations fo- 
cus on minimizing the magnetic energy in specific, chosen 
sublattice structures, other possible commensurate or in- 
commensurate states are effectively neglected. Calcula- 
tions on larger magnetic unit cells or triangular lattices 
with periodic boundary conditions may uncover larger 
sublattice or more complex spin configurations within 
the model that may be relevant, as was found to be the 
case for the zero-field phase of doped CuFe02.'^^'^ The 
recently proposed incommensurate 120°-like spin struc- 
ture above sjlg is one example, though its underlying 
Landau theory does not capture some general features 
of the exp erimental high-field magnetization curve at 
present .123 A fully accurate and quantitative description 
of the magnetism of CuFe02 would require the inclusion 
of all additional features of the system that could play a 
role. The incorporation of finite temperature, a more re- 
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FIG. 11: (Color online) Calculated minimum energy per spin for the commensurate sublattice phases [B _L c) for = Q (a)) 
and Jz = 0.25Ji (b)). Dashed vertical lines indicate experimental transition fields. Inset: Corresponding field dependence of 
the simulated magnetization (blue) and assumed easy axis anisotropy (red) in CuFe02 for _B _L c. 



alistic pho non m odel (yielding longer range biquadratic 
interaction^^mi) and quantum spin effects may improve 
the quantitative understanding of the system. Further- 
more, more exotic interactions may play a role in stabiliz- 
ing the incommensurate spiral state."^ Nevertheless, the 
simple PCS model presented here is shown to capture 
almost all general features of the experimental magne- 
tization process in both field configurations, providing 
a satisfactory and intuitive description of the observed 
magnetism in CuFe02. 

IV. CONCLUSIONS 

Summarizing, we have performed magnetization 
experiments on CuFe02 at various temperatures below 
Tn2 up to high magnetic fields, both for B \\ c and 
i? _L c field configurations. The characteristic magnetic 
staircase of CuFe02 was reproduced and found to retain 
its general features with increasing temperature below 
T/V2- As the temperature approaches Tn2 however, 
transition features are progressively smoothed out and 
plateau phases are found to acquire increasing slopes. 
Moreover, the transition from the coUinear 5SL to the 
coUinear SSL" phase (at -BI3) was shown to split up into 
a two-step transition near T/V27 revealing an additional, 
possibly noncollinear, intermediate state / at these 
temperatures. Additionally, the various critical fields 
of the same nature are shown to exhibit a very similar 
temperature dependence; all first order transitions 
exhibit an analogous relative decrease with temperature, 
and second order transitions are found to be relatively 



temperature independent. Correspondingly, we have 
thoroughly mapped out the experimental _B,T phase 
diagrams of CuFe02 for both the parallel [B \\ c) and 
perpendicular {B _L c) configurations and expanded 
them in both temperature and magnetic field. Through 
numerical minimization of the PCS model applied to the 
consecutive commensurate sublattice phases of CuFe02, 
also the low-field part of the experimental magnetization 
process was adequately simulated, yielding reasonable 
estimates for the additional parameters J2 and J3. 
Incorporation of an additional interlayer exchange 
interaction in the model was shown to result in a nearly 
identical simulation and a somewhat adapted set of 
exchange interactions. Thus, the proposed PCS model, 
combined with the underlying notion of progressive 
symmetry increase with applied field, is found to provide 
a satisfactory semi-quantitative description of the entire 
magnetization process of CuFe02. 
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